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^ ' melting crystal picture and the amplitudes are given by the quantum dilogarithm. In the 
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show that these two amplitudes are related by a certain integral transformation. We argue 
that this transformation represents the deformation of closed string background due to the 
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1. Introduction 

The gauge/string duality is one of tlie most important aspects of string theory. It 
appears that this duality is also essential for the topological string theory. In [mi, it is 
realized that the 't Hooft expansion of Chern-Simons theory on is exactly the same as 
the topological A- mo del on the resolved conifold 1)©(9(— 1)-^P^. It was shown that 
the partition function of SU (N) level k Chern-Simons theory on ^ 

1 -i-r vrfa, n) 

Zcs = - TTn^n^smff^ 1-1 

agrees with the partition function of topological closed string on local [ffl 



oo 



- n ( ^ (1-2) 

n=l ^ ^ / 

up to a non-perturbative factor r]{q)^ . In this correspondence, various parameters are 
related as 

k + N ' ' (1.3) 

q = e-3^ , Q = e-\ 

where Qs is the string coupling and the 't Hooft parameter t is identified as the Kahler 
parameter of P^. This is an example of geometric transition where D-branes wrapped 
around in the deformed conifold [H is replaced by a S-field on S"^ in the resolved 
conifold. In this duality was derived by using a linear sigma model. 

In this paper, we consider open string amplitudes associated with non-compact D- 
branes in this background. There are two ways to define this amplitude. The first way 
is to introduce non-compact Lagrangian A-branes in the deformed conifold T*S^. The 
amplitude of the A-branes is given by Wilson loops in Chern-Simons theory 0. Later it 
was shown that the same amplitude is obtained in the melting crystal picture . In this 
picture, the A-brane corresponds to a defect of crystal and it is represented by a certain 
fermion operator "^foix) [|T0| . The amplitude {"i! d{x)) is given by the quantum dilogarithm 

EMI- 

The second way to define the D-brane amplitude is to use the matrix model description 
of Chern-Simons theory p^^T7|. A natural D-brane amplitude in the Chern-Simons matrix 



model is the expectation value of determinant (det(a; — M)). Since this matrix model 



appears as the mirror B-model of conifold [|T5[, the determinant det(x — M) represents a 



1 



non-compact B-brane in the mirror description |T^. It is well known that (det(a; — M)) 
is equal to the A^-th orthogonal polynomial Pn{x). For the Chern-Simons matrix model, 
the associated orthogonal polynomial is known as the Stieltjes-Wigert polynomial ||19|| . 

It is natural to ask what is the relation between the A-brane amplitude in 
the crystal picture and the B-brane amplitude (det(a; — M)) in the Chern-Simons matrix 
model. In this paper, we will show that they are related by a simple Gaussian integral. 

This paper is organized as follows. In section 2, we review the Gaussian matrix 
model as a preliminary of the analysis of Chern-Simons matrix model. In section 3, we 
find a simple relation between A-brane amplitudes and B-brane amplitudes. Section 4 is 
discussion. 



2. Gaussian Matrix Model: Review 

Before discussing the D-branes in topological string, let us first review a simple ex- 



ample, namely FZZT-branes in Gaussian matrix model [^. Although this is a too trivial 
example, we will see that most of the properties found in the Gaussian matrix model have 
close analogues in the Chern-Simons matrix model. 
The Gaussian matrix model is defined by 



Zg= I dMe-^^^\ (2.1) 

J NxN 

We are interested in the expectation value of the determinant operator det(a; — M) 



1 



t2 



(det(x - M)) = ^ / dMe~^^'^ det{x - M), (2.2) 
Zg 



which can be interpreted as the wavefunction of FZZT-brane. As discussed in detail in [|2C 



this integral ( p.2|) is evaluated by rewriting the determinant as a fermion integral and then 
integrate out M. The resulting expression is an effective theory on the single FZZT-brane 

ds 12 

(det(a: - M)) = / —=6"^^ fN{x + is). (2.3) 



Here we introduced the function /Ar(x) = . We can easily see that (|2.3|) is nothing but 
the integral representation of the Hermite polynomial, which in turn is the A^-th orthogonal 
polynomial Pn{x) of Gaussian measure 

(det(x - M)) = P^{x) = (I) ' Hr, (^-^ j . (2.4) 



It is interesting that fN{x) = can be thought of as a "classical" value of the 
determinant. Namely, fN^x) is the value of det(x — M) evaluated at the minimum of the 
Gaussian potential TrM^ 

/7v(x) = det(x-Mo) , Mo = diag(0,0,---,0). (2.5) 



Prom this viewpoint, ( |2.3| ) defines an integral transformation Q which maps the "classical 
wavefunction" fN^x) into the "quantum wavefunction" Pn{x) 

Q : det(x - Mq) (det(a; - M)) 

Note that this relation holds for arbitrary N , therefore fn{x) = x"^ is mapped to Pn{x) for 
all n 

Q:U{x)^Pn{x). (2.7) 

We can generalize this picture to the multi-point correlators of FZZT-branes. Follow- 
ing the similar procedure as above, the ii'-point function of determinants is written as a 
K X K matrix model E 



rrdet(x, -M)\ = /" dS e-^"^'^' det{X + iS)^, (2.8) 
u=i / -^^x^ 

where X = diag(xi, ■ ■ ■ , xk)- This can be thought of as an effective open string theory 
on the K FZZT-branes.0 After integrating out the angular part of 5", (2TB) becomes an 
integral over the K eigenvalues of S 

A(x) (jldetix, -M)^= I II -^e-^^'/W(xi + is^,...,XK + isk). (2.9) 
Here A{x) = Y[i>ji^i ~ ^j) the Vandermonde determinant and f^^\x) is given by 

K 

(xi, ■ • • , xk) = A(x) n fN{x,). (2.10) 

i=l 



^ After taking a double scaling limit, this open string theory on the FZZT-branes is identified 



as the Kontsevich model [21,20|. 



As before, we regard f^^\x) as a "classical" value of the X-point function. The true 
correlator of determinants is given by the Q-transform of f^^\x) (^). 
To evaluate the integral ( |2.9| ), we rewrite f^^\x) ( p.lO| ) as 



K 



= det(a;r') = det(a;f+^-^) = det(/^+,_i(x,)). (2.11) 



2=1 



Now the integral ( |2.9| ) is easily evaluated by recalling that the integral transformation Q 
maps fn{x) into Pnix) ( p.7| ). Therefore, the i^-point correlator of FZZT-brane is given by 



K 



A{x) ( ndet(a:, - M) ) = det(P^+,_i(x,)). (2.12) 



\i=l 



As expected, this agrees with the general formula of the correlator of determinants 
In ( ^.12| ) we have multiplied the Vandermonde determinant A(x) to the correlator of 
determinants. This factor makes the left hand side of ( ^.12| ) anti-symmetric in Xi, therefore 
FZZT-branes become fermionic To summarize, we found the following relation 



K IKK 

Q: A(x) JJdet(x, -Mo) ^ A(x) / JJ det(xi - M) \ . (2.13) 

i=l \i=l I 

In the next section, we will see that the transformation Q is exactly what we are 
looking for, i.e., Q maps A-brane amplitudes into B-brane amplitudes. 



3. D-brane Amplitudes on Conifold 

3.1. A-branes in the Calahi-Yau Crystal Picture 

As shown in 1^,0, the target space theory of topological A-model is reformulated as 
a certain U{1) gauge theory and it leads to the quantum foam picture of Kahler gravity, 
which in turn is related to the statistical model of crystal melting. In the case of resolved 



conifold, the relevant crystal is identified in |Tl|]. It is given by a lattice in positive octant of 
M'^ with lattice spacing and it is restricted to the interval [0, Nqs] in one direction. This 
statistical problem is solved by slicing the 3D partition into a sequence of 2D partitions 
and use the relation between 2D Young diagram and free fermion 0]. 

In this crystal picture, a non-compact A-brane is represented by a half-line in R'^ and 
it creates a defect in the lattice . It is further shown that the presence of defect is rep- 
resented by the insertion of a certain fermion operator \E' in the free fermion computation 



4 



mentioned above. If the half hne associated with the i-th D-brane ends on one of the axis 
of Ri^ at tti = gs{Ni + |), we identify the moduh Xi of this D-brane as i 

x. = g^^ (3.1) 

We are interested in the amphtude of Lagrangian A-brane with moduh x 

Zn{x) = {^d{x)). (3.2) 

This is first obtained by a Wilson loop calculation in Chern-Simons theory on S""^ 0, and 
later it is reproduced by the melting crystal picture [lT0| . [Tl| , |T3| . We will consider an A-brane 
with the corresponding half-line ending on the interval [0,Ngs]. Then the A-brane has a 
topology X [p3tp4| , |To[| . The explicit form of the single brane amplitude is given by 

where L{x, q) is the quantum dilogarithm ^5 



Lix,q)=ll{l-xqn. (3.4) 

n=l 

Similarly, the amplitude in the presence of K A-branes is given by |[T0| , |TTl , p!3[] 

K K K 

Z{xi, ...,xk) = X{{1- x,x-^) J] Zi,{xi) = A(a;) J] Z^(xO J] x]-\ (3.5) 

i<j i=l -1=1 i=l 

Although it is a little bit ad hoc, it seems natural to drop the last factor Y\.i x\~^ in order to 
make the multi-point function of D-branes anti-symmetric under the interchange of XiM. 
This is motivated by the intuition that non-compact D-branes are fermions. Then the 
X-point function of D-brane operators becomes 

K \ K 

\[^D{Xi) \ = ^{x)\[ZN{Xi). (3.6) 

\i=l I i=l 

One can immediately notice the similarity of this expression to the "classical" correlation 
function ( p.lO| ) in the Gaussian matrix model. We would like to argue that this is not 
just a coincidence. We will see below that the "quantum" version of (|3.6|) is given by the 
correlator of determinants in the Chern-Simons matrix model. 



9 — 

This definition of Xi defers from the usual convention by a factor of g 2 . Our definition is 
convenient when comparing to the B-brane amplitude in the next section. 

^ Perhaps this factor might be canceled by carefully analyzing the zero-mode of boson appear- 
ing in the bosonization of fermion. 
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3.2. B-branes in Chern-Simons Matrix Model 

In this subsection, we will consider the D-brane amplitudes in the Chern-Simons 



matrix model. In [14|, it is observed that the partition function of Chern-Simons theory 



is written as a matrix integral. It is then realized that this Chern-Simons matrix model 



naturally appears in the B-model mirror of resolved conifold [|T^ . 

Let us recall the derivation of Chern-Simons matrix model. Using the Weyl denomi- 
nator formula, it is easy to see that the partition function of Chern-Simons theory ( |1 . 1| ) is 
written as an integral 



Zcs= Yldu^lli 2 smh^^^^] e-^^y^. (3.7) 
i=i i<j ^ ^ 

This is almost a matrix model integral, but the measure factor is not the usual Vander- 



monde determinant. By the following change of variables |2f:,16 



m, = e'^'q-'^ (3.8) 

we can change the measure factor in (|3.7|) into the usual Vandermonde determinant A(m)^. 
Then Zcs becomes an x hermitian matrix model with log-squared potential 

Zcs = j dMe-^^'^'°^^'^\ (3.9) 



Now let us consider the D-branes in this model. From the general argument [|T| 
non-compact B-branes correspond to the insertion of determinant operators in the matrix 
model. We define the B-brane operator with moduli x as 

B{x) = {-l)^q^'+^^ det{x - M) = det{q^+^M - q^+^x). (3.10) 

We have put a g-dependent prefactor for later convenience. This is partly motivated by the 
relation (|3.8[) that log(g^mi) is the natural Gaussian variable Ui in the original integral 
( ^77|) . The additional factor q^ is related to our definition of D-brane moduli x (see footnote 
2). Then the B-brane amplitude is given by 

Sn{x) = {B{x)), (3.11) 



where the expectation value is taken in the Chern-Simons matrix model ( |3.9| ). From the 
Heine's formula, the expectation value of determinant is the A^-th orthogonal polynomial. 
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In the case of log- normal measure ( p.9|) , the orthogonal polynomial is known as the Stieltjes- 
Wigert polynomial [[T^B With our normalization of B{x) ( ^.101 ), the A^-th Stieltjes-Wigert 
polynomial is given by 

AT _ 

N 
k 



Sn(x) 



N 



A;=0 



(3.12) 



Here 



N 
k 



denotes the g-binomial which is defined by 

i<l)N 



N 
k 



(?)n=n(i-^')- 



(3.13) 



{q)k{q)N-k ' 

We can also consider the multi-point function of B-branes. From the general formula of 



correlator of determinants , the ii'-point function of B-branes reads 

K 

det{SN+j-i{xi)). 



(3.14) 



\i=l 



3.3. A Map from A-brane Amplitudes to B-brane Amplitudes 



Now we consider the relation between the A-brane amplitude Zn{x) = (\E'd(x)) in the 
crystal picture and the mirror B-brane amplitude Sn{x) = {B{x)) in the Chern-Simons 
matrix model. To see this relation, let us rewrite Z]\[{x) given in the product form 
into a summation form ||19|| i 



N r 



k=0 



N 
k 



,ifc(fc+i) 



-X] 



(3.15) 



By comparing this expression with Sisf{x) in ( |3.12| ), one can immediately see that they 
are almost identical except for the exponent of q. This difference is taken care of by the 
following integral 

/7 Q 19.. 1 . O 

(3.16) 



Therefore, Sn{x) and Zn{x) are related by 



1_„2 

-e ^3s e 



q^" . 



-.e~^''ZN{xe''). 



(3.17) 



In [11 1, a unitary matrix model with measure i?3(g, e**^) is considered. The orthogonal poly- 
nomials on 5^ associated with this weight are known as the Rogers-Szego polynomials |19]. They 
are related to the Stieltjes-Wigert polynomials by a change of variable. 
^ This is a special case of the g- binomial theorem |27]. 
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As advertised, this is nothing but the integral transformation Q appeared in the Gaussian 
matrix model! (Recall that the variable obeying the Gaussian distribution is log a; in the 
Chern-Simons matrix model. Therefore, the transformation (|3.17|) is exactly the same as 
Q in ( |2.6|) .) To summarize, the A-brane amplitude and the B-brane amplitude are related 
by the map Q 

Q:Z^{x)^S^{x). (3.18) 

We can easily generalize this relation to the multi-point correlators. Let us consider 
the 2-point function for illustration. The 2-point function of A-branes is (see ( ^.6| )) 

Y[ <^fD{x^y) = {X2 - Xi)Zn{xi)Zn{x2). (3.19) 

Using the recursion relation of ^^^(x) 

xq^'+^ZNix) = Zn{x) - Zn+i{x), (3.20) 

( p.l9|) is rewritten as 

n ^d{x{) \ = q-^-^ [ZN+l{xi)ZNix2) - ZNixi)ZN+l{x2)] . (3.21) 
\i=l I 

By applying the Q-map ( |3.18| ) , this becomes the 2-point function of B-branes ( |3.14|) , up to 
an overall factor q~^~^ . One can easily generalize this argument for the general i^-point 
function, and it is found to be 

JJ ^D{xi) \ = det{ZN+j-i{x,)). (3.22) 
u=i / 

Here we suppressed an overall q dependent factor for simplicity. Clearly, (|3.22| ) is mapped 
to ( |3.14| ) by the Q-map (|3.18| ). In this way we find that the A-brane correlator and the 
B-brane correlator are related by the Q-map 

A(x)^ni?(x.)^= y ^n72^^""^'^\^n^^(^^-^'^^)/- (3-23) 

Recalling the form of A-brane correlator ( |3.6| ), this relation is written in a similar form as 
the Gaussian matrix model case ( |2.13D 

Q: /\{x)\[ZN{xi) ^ /\{x) {\[B{xi)\ . (3.24) 

1=1 \i=l / 



Finally, if we introduce the Q-transform Q'i'D of A-brane operator ( |3.23| ) can be 
written in a more suggestive form 



K \ IK 



A(x) ( n B{xi) ) = ( n 2*^(^0 ) ■ (3.25) 



\i=l I \i=l 



Note that d{x) is smeared by the Gaussian integral with width Qs, thus it is no longer 
a local operator in the x-space. Perhaps it might be related to the noncommutativity in 
spacetime. 



4. Discussion 

In this paper, we found that the A-brane correlators and the B-brane correlators are 
related by the transformation Q, which is just a Gaussian integral. We would like to 
understand the physical meaning of this relation. By the analogy to the Gaussian matrix 
model case, it seems natural to regard the A-brane amplitude Zjq{ QiS Oil classical value 
of the determinant in the Chern-Simons matrix model. From the product form of Zn{x) 
(p73|), the classical background Mq corresponding to Z]\[{x) is given by 

Z^ix) = (-l)^gi^(^+i) det(x - Mo), 
Mo = diag((7-\Q-2,---,(7-^). 

In other words, Zn^x) vanishes at 

logx = gs,2gs,---,Ngs, (4.2) 

i.e., the zeros of Zn{x) are equally spaced.! This is consistent with the picture that the 
A-brane probes the background lattice of the crystal. Namely, the A-brane amplitude 
obtained in the crystal picture is a probe approximation ignoring the backreaction. More- 
over, it is known [18,1^ that the A-brane amplitude is the zero energy wavefunction of the 



Hamiltonian H{u, —ggdu), where H{u, v) defines the mirror Riemann surface H{u, v) = 0. 
In our case, ZNie"^) is the zero-mode of the Hamiltonian 



H{u,v) = l-e'' -{l-e''q^)e\ (4.3) 



^ In [28 1, the eigenvalue density of Chern-Simons matrix model is studied numerically at finite 
A''. It is observed that the density has peaks with almost equal spacing. This is consistent with 
our identification of Mq as the classical configuration of Chern-Simons matrix model. 
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This agrees with the known mirror of the resolved conifold [15,2^ up to a Unear canonical 
transformation of the u,v coordinates ||18|| . 

On the other hand, the B-brane amplitude Sn{x) includes the effect of deformation 
of closed string background due to the presence of D-brane Jl^. In fact, the insertion of 
determinants is equivalent to a certain shift of matrix model potential |lT8|,pO 



JJdet(xi - M) 



1 



Zcs 



dM exp 



2g. 



— Tr(log M)2 + V Tr log(x. 



M) 



(4.4) 



As argued in [18|, the deformation of matrix model potential due to the insertion of B- 
branes corresponds to turning on gravitational descendants of Kahler class of P^. It is 
remarkable that this deformation is captured by the simple map Q. Another way to see 
this effect is to look at the zeros of Sn{x). For example, S2{x) is factorized as 



<S2ix) = 1 — g2 (1 -f q'jx + q X = q {x — a_|_)(x — a 
1 



a± = -g 2 



l + q±^{l-q)il + 3q) 



(4.5) 



The zeros are no longer equally spaced. This might be seen as the distortion of crystal due 
to the insertion of brane. One might expect that the "quantum" wavefunction Sn{x) is 
significantly different from its classical counterpart Zn{x). It would be interesting to study 
the large behavior of Sn{x). It would be also interesting to consider a double-scaling 
limit of Chern-Simons matrix model as in the case of Gaussian matrix model . 

In this paper we did not consider the insertion of anti-brane. In the crystal picture, 
it is known how to construct the anti-brane operator However it is not clear to 

us what is the corresponding operator in the Chern-Simons matrix model. We leave this 
as a future problem. 
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